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ABSTRACT 


# 

Turbulent  end  wall  boundary  layers  in  the  hub  region  of  axial,  com¬ 
pressor  blade  rows  have  been  studied  in  the  Gas  Turbine  Laboratory. 

* 

Both  analytical  and  experimental  work  are  presented.  An  attempt  was 
made  to  predict  the  growth  of  the  boundary  layer  momentum  thickness  as 
it  passes  through  a  cascade  of  blades  using  the  momentum  integral 
equation  and  several  assumptions.  This  calculation  was  compared  with 
measurements  taken  along  an  assumed  free  streamline  in  a  stationary 
cascade  of  blades  supplied  with  an  artifically  skewed  inlet  boundary 
layer.  From  this  comparison,  the  relative  influence  of  individual 
terms  in  the  momentum  integral  equation  is  deduced.  The  effects  of 
inlet  skewing  and  cross  flow  in  the  boundary  layer  in  the  region  between 
blades  is  described.  Two  methods  of  representing  main  and  cross  flow 
velocity  profiles  sure  compared  with  experimental  data. 
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Introduction 


The  flow  near  the  casing  of  axial  compressors  is  difficult  to  pre¬ 
dict  since  viscous  forces  are  appreciable.  If  the  flow  near  the  wall 
at  the  exit  of  a  row  of  blades  could  be  calculated,  given  the  inlet 
flow  conditions,  the  blade  geometry  and  the  end  wall  configuration, 
single-stage  and  multi-stage  compressor  performance  could  be  more  re¬ 
liably  estimated.  The  problem  can  be  considered  as  one  of  the  boundary 
layer  type,  the  turbulent  case  being  most  important  for  compressors . 

A  study  directed  toward  the  solution  of  the  problem  has  been  conducted 
in  the  M.I.T.  Gas  Turbine  Laboratory.  Some  interesting  results,  mostly 
of  a  qualitative  nature,  have  come  from  the  work  and  are  reported  in 
this  paper. 

2  Description  of  Flow  in  the  Hub  Region 

The  blade  configuration  at  the  casing;  i.e.  clearance  between  the 
tips  of  the  blades  and  the  wall  and  any  motion  of  the  casing  relative 
to  the  blades,  strongly  influences  the  character  of  the  problem.  Since 
the  inlet  flow  to  a  given  row  is  determined  by  the  previous  row,  the 
configurations  of  the  boundaries  in  two  successive  rows  must  be  speci¬ 
fied  for  any  particular  problem.  We  have  limited  our  investigation  to 
one  overall  configuration., of  an  axial  compressor  blade  row  where: 

1.  No  clearance  exists  between  the  blade  tip  and  the  casing. 

2.  The  easing,  or  end  wall,  is  Stationary  with  respect  to 
the  blades . 

3«  The  previous  row  has  the  same  configuration. 

This  overall  configuration  would  be  found  at  the  hub  of  a  compressor 
with  shrouded  stator  roots  as  shown  in  Figure  1. 

Consider  the  flow  at  Section  A-A,  Figure  1*  The  velocity  triangles 
at  that  point  might  be  as  shown  below. 
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Velocity  Profiles  in 
x'  and  z1  Directions 
Relative  to  the  Rotor 


Velocity  Profiles  in 
x  and  z  Directions 
Relative  to  the  Rotor 


Primes  denote  quantities  relative  to  rotor. 


By  addition  of  the  wheel  speed  vector  (rf X  )  to  flow  relative 
to  the  stator  one  gets  the  flow  relative  to  the  rotor.  It  is  inter¬ 
esting  to  notice  how  a  boundary  layer  which  is  essentially  unidirectional 
(such  as  the  one  shown  leaving  the  stator  above)  can  become  highly 
skewed  relative  to  the  rotor  by  adding  the  wheel  speed  vector.  In  the 
above  figure,  U  refers  to  the  exit  velocity  from  the  stator  and  U*  to 
the  entrance  velocity  to  the  rotor,  both  at  some  distance  away  from  the 
casing,  while  V  and  V'  denote  the  same  velocities  at  seme  point  in  the 


3. 

boundary  layer  near  the  inner  casing.  The  flow  entering  the  rotor  blades 
row  near  the  wall  (i.e.  the  entering  boundary  layer)  has  a  high  angle 
of  incidence  on  the  blades  (which,  we  assume,  have  been  designed  to 
accept  the  free  stream  angle)  and,  unlike  conventional  boundary  layers, 
has  kinetic  energy  comparable  to  the  free  stream.  The  velocity  varies 
from  r_TL  ,  the  wheel  speed  directed  tangentially,  at  the  rotor  hub 
to  U',  directed  at  some  design  angle  (in  the  neighborhood  of  4o°  or  50°, 
measured  from  the  axial  direction)  at  sane  distance  away  from  the  wall. 

As  the  boundary  layer  fluid  flows  on  to  the  rotor,  high  shear 
stresses  near  the  wall  (Figure  l)  dissipate  some  of  its  kinetic  energy. 

The  boundary  layer  fluid  progresses  through  the  blade  row,  acted  upon 
by  pressure  forces  add  viscous  forces,  with  changing  velocity  and 
direction.  The  pressure  field  tends  to  accelerate  any  velocity  de¬ 
ficient  fluid  toward  the  suction  side  of  the  blade  passage  and  a  flow 
as  shown  in  Figure  2  is  common. 

Thfe.  same  events  occur  as  the  fluid  progresses  through  the  succeeding 
stator . 

5  The  Investigation 

3.1  General  Approach  and  Some  Assumptions 

The  discussion  will  be  limited  to  an  incompressible  fluid.  We  con¬ 
sider  the  flow  relative  to  the  blades' and  will  refer  to  the  flow  out¬ 
side  the  boundary  layer  as  the  free  stream.  The  change  in  the  radius, 
measured  from  the  axis  of  the  machine  to  a  free  streamline,  will  be 
assumed  small  so. that  the  stagnation  pressure  based  on  the  relative 
velocity,  p  +  p|ya .will  be  constant  along  a  free  streamline.  The  region 
2 

of  shear  flow  near  the  casing  will  be  considered  a  proper  boundary  layer. , 
with  small  thickness,  which  is  turbulent. 

It  is  convenient  to  think  in  terms  of  a  curvilinear  coordinate 
system  where  the  x-axis  is  in  the  direction  of  the  free  stream  just  out¬ 
side  the  boundary  layer,  the  y-axis  is  perpendicular  to  the  end  wall  or 
casing  and  the  z-axis  is  perpendicular  to  the  x  and  y  axes  so  as  to 
form  a  right  hand,-  orthogonal  system.  The  velocity  of  the  free  stream, 
just  outside  the  boundary  layer  is  U  and  the  velocity  components  in 
the  boundary  layer  along  the  x,  y  and  z  axes  are  u,  v  and  w  respectively, 
w  is  often  called  the  "cross  flow"  velocity. 

The  blades  will  determine  the  flow  field  just  outside  the  boundary 
layer  and  since  we  are  primarily  concerned  with  the  boundary  layer  flow, 
we  assume  the  velocity,  static  pressure  and  curvature  at  each  point  of 
every  free  streamline  are  known.  For  thin  boundary  layers,  £j?  q 

ay  ~ 

and  the  pressure  through  the  boundary  layer  is  determined  by  the  free  stream. 
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A  complete  solution  of  the  problem  would  yield  the  velocity  and 
direction  of  flow  within  the  boundary  layer  over  a  region  between  two 
blades  such  as  the  one  in  Figure. 2c  This  information  could  also  be 
presented  in  terms  of  the  boundary  layer  flow  along  several  different 
x-axes,  each  being  a  different  free  streamline.  In  our  study,,  attention 
has  been  focused  on  just  one  free  streamline .axis,  the  one  which  leaves 
the  blade  row  approximately  halfway  between  two  blades ,  Figure  2c 


3 . 2  Some  Observations  in  a  Cascade  of  Compressor  Blades 


5.2.1  Simulation  of  Flow  in  Compressor  Blade  Row  in  a 
’’"Rectilinear  Cascade  ~~ 


Studying  the  behavior  of  the  wall  boundary  layer  in  an  actual 
machine  is  quite  difficult  for  several  reasons .  Measurements  relative 
to  the  blades  at  inlet,  between  the  blades  and  at  outlet  are  necessary. 
In  a  rotor  such  measurements  involve  considerable  mechanical  com¬ 
plication.  In  a  stator  following  a  rotor,  fluctuations  of  the  flow 
with  time  are  present.  In  addition,  the  dimensions  of  the  flow  passage 
between  blades  in  a  machine  are  generally  too  small  to  permit  accurate 
measurements . 


To  avoid  these  difficulties,  a  cascade  was  constructed*  to 
simulate  the  flow  near  the  hub  of  an  axial  compressor  rotor.  An  axial 
entry  rotor  was  simulated  to  simplify  future  comparisons  with  measure¬ 
ments  in  an  actual  machine.  (Some  measurements  have  been  made  in  an 
axial  entry- rotor  but  no  comparisons  have  yet  been  possible.) 


The  axial  entry  rotor  meets  the  conditions  described  in  Section  2. 
and  contains  the  important  features  of  hub  flow,  the  dkewness  and 
relatively  high  kinetic  energy  of  the  boundary  layer  at  the  blade  row 
entrance . 


The  testing  of  a  cascade  with  skewed,  high  kinetic  energy  inlet 
flow  near  the  wall  poses  several  new  problems  so  the  initial  investiga¬ 
tion  was  exploratory  and  made  use  of  available  equipment.  As  a  result, 
the  accuracy  of  measurements  (particularly  flow  direction)  was  not 
very  good.  The  apparatus  and  measuring  techniques  are  described  in 
Appendix  A. 


3.2.2  Blade  Geometry,  Free  Stream  Properties  and  Inlet  Conditions 


The  airfoils  had: 


Circular  arc  camber  line 

NACA  Four  Digit  Series  thickness  distribution  with  9$ 
maximum  thickness 


*  The  construction  and  testing  of  the  cascade  were  the  subjects  of  an 
S.M.  Thesis  in  the  Gas  Turbine  Laboratory,  Reference  5* 
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i.' 
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Chord j,  C  =  2.8" 

Camber  angle,  9  =  35° 

The  cascade  was  set  within: 

Stagger  angle,  =34.3° 

Pitch-chord  ratio,  s/c  =  0.8g8 

The  free  stream  flow  entered  at  =  52°  and  left  at  p2  =  22.5° 
(angles  measured  as  shown  in  Figure  3)»  The  resulting  distribution  of 
static  pressure  at  the  wall  is  shown  in  Figure  4.  The  distribution  of 
dynamic  head,  q  (q  =  p  -  p),  along  the  axis  shown  in  Figure  4  (assumed 
to  be  a  free  streamline;  is  plotted  in  Figure  5°  The  points  B,  E,  G, 

H,  I  and  J2  refer  to  measuring  stations  along  the  center  streamline  as 
shown  in  Figure  6.  The  overall  static  pressure  rise  coefficient,  Cp, 
for  the  free  stream  was  Cp  =  0.337=* 

The  inlet  conditions  to  the  axial  entry  rotor  were  measured  in 
the  machine,  then  approximated  in  the  cascade.  The  inlet  velocity,  V, 
relative  to  the  stationary  hub  is  plotted  versus  y,  the  distance  from 
the  hub,  in  Figure  7=  It  should  be  noted  that  in  this  case  the 
boundary  layer  relative  to  the  stationary  hub  was  not  skewed  and  in 
the  axial  direction.  To  get  the  flow  in  the  streamline  coordinate 
system  relative  to  the  moving  blades,  a  vector  addition  of  the  wheel 
speed  to  the  air  velocity  relative  to  the  casing  was  performed  at  each 
value  of  y.  The  relative  velocity  was  reduced  to  its  components  u  and 
w,  as  shown  in  the  diagram. 


Free  streamline 
direction  relative 
to  rotor  blades 


*  This  value  differs  from' the  value  predicted  by  two-dimensional  theory 
because  the  cascade  did  not  behave  uniformly.  One  of  the  five  blade 
passages  was  badly  stalled  over  its  entire  span. 
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where:  U  is  velocity  relative  to  the  casing  outside  the  boundary  layer 

V  is  velocity  relative  to  the  casing  at  a  point  within  boundary  layer 

U'  is  the  velocity  relative  to  the  moving  blades  at  y  =  5 

V!  is  the  velocity  relative  to  the  moving  blades  within  the  boundary  layer 

x  and  z  are  streamline  coordinates  relative  to  the  blades 
u  and  w  are  velocity  components  in  the  x  and  z  directions 
=  1.075  for  the  case  simulated. 

Y'.-TL 

The  resulting  inlet  velocity  profiles  are  plotted  in  Figure  8 
along  with  the  cascade  approximations . 

5.2.3  Boundary  Layer  Flow  Through  Blades  and  Exit  Conditions 

The  u  and  w  velocity  profiles  were  measured*  at  points  A,  B,  C, 

D,  E,  F,  G,  H,  I,  Jx,  J2,  J3  along  the  bottom  wall  of  the  cascade, 

Figure  6,  and  the  results  are  presented  in  Figures  10  and  11.  The  progress 
of  the  boundary  layer  flow  along  the  free  streamline  axis  can  be  ob¬ 
served  from  the  measurements  at  B,  E,  G,  H,  I  and  J2.  Some  idea  of  flow 
variations  across  the  passage  can  be  obtained  from  the  other  profiles. 

The  integral  quantities  appearing  in  the  momentum  integral  equations  ** 
have  been  calculated  from  the  velocity  profiles  and  are  tabulated  in 
Table  I. 


The  bottom  wall  of  the  cascade  was  coated  with  a  mixture  of 
kerosene  and  carbon  black  before  one  of  the  runs .  When  the  flow  was 
turned  on,  traces  were  made  on  the  coating  which  are  roughly  indicative 
of  the  flow  direction  adjacent  to  the  wall.  Figure  12  is  a  photograph 
of  the  resulting  traces.. 

5.2.4  Blade  Forces 


The  static  pressure  distribution  around  the  blades  was  measured 
at  several  distances,  Y,  from  the  end  wall.  The  axial  and  tangential 
components  of  the  pressure  force  on  the  blade  were  computed  and  are 
plotted  vs.  Y  in  Figure  15 .  The  forces  are  normalized  with  the  value 
at  Y  =  1.5  inches,  which  was  taken  to  be  the  free  stream  condition. 

5.5  Calculation  of  9x  Along  Assumed  Streamline  Axis  *** 

The  momentum  integral  equation  along  a  free  streamline  axis  is: 


as* 


+  iLtir)  4. 

fr  a  / 


dz 


(1) 


*  Measuring  techniques  are  described  in  Appendix  A.  ~ 

**  See  Appendix  B  for  momentum  integral  equations . 

***  For  a  discussion  of  calculation  methods,  see  Appendix  B. 
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Mager  (Reference  l)  was  able  to  calculate  the  growth  of  along  a 
free  streamline  axis  in  a  90°  bend  by  neglecting  the  term  contributed 
by  the  presence  of  cross  flow,  •  Using  this  assumption  for 

b7- 

the  present  case,  the  variation  of  ©x  along  the  assumed  streamline  was 
calculated  with; 


1. 

2. 

3» 


The  measured  variation  of  q  with  x,  Figure  5* 

H  =  1.30,  an  average  value  suggested  by  the  data 
The  Squire  and  Young  formula  for  skin  friction  as  presented 
in  Reference  3>  ^  n ^ 


£*  *  HJ-*-07s  *•) 

Lox.  *-  lo 


4.  The  initial  from  the  data  at  point  E,  9zg  =  0 .0437  inches 


The  calculations  were  started  at  point  E  and  carried  to  point  J2> 
no  attempt  was  made  to  include  the  high  shear  region  from  B  to  E.  The 
result  of  the  calculation  is  plotted  in  Figure  l4. 


4  Discussion  of  Results 

4.1  Qualitative  Description  of  the  Boundary  Layer  Flow  Along 
the  Assigned  Streamline  Axis 

At  the  leading  edge  of  the  splitter • plate  the  velocity  is  high  and 
is  directed  parallel  to  the  plane  of  the  cascade.  (Simulating  the  flow 
on  to  the  rotor  hub  where  the  velocity  adjacent  to  the  rotor  hub  is 
Y'YJl directed  tangentially.)  This  velocity  is  immediately  reduced  to  zero 
and  the  velocities  near  the  wall  quickly  decreased  in  a  region  of 
high  shear  stress  from  B  to  ]£§!  Figure  6.  This  shearing  action,  combined 
with  the  decelerations  caused  by  the  pressure  field,  causes  a  decrease 
in  the  z  component  of  velocity,  w,  and  sin  increase  in  the  u  component 
near  the  wall  from  B  to  E.  As  a  result,  the  angle  between  the  stream¬ 
line  at  the  wall  and  the  free  stream  goes  from  24.5°  at  B  to  10.0°  at  E. 
This  pronounced  turning  near  the  wall  is  indicated  by  the  carbon  black 
traces,  Figure  12.  It  is  interesting  to  note  that  in  this  region  ©x 
is  decreasing  in  the  x  direction,  though  the  x  pressure  gradientsis 
adverse  (pressure  rising  in  the  x  direction). 

The  component  of  the  pressure  gradient  in  the  x  direction, 

is  positive,  while  the  gradient  in  the^d:id4pjfedtl9n , ,  is  negative 

at  all  points  from  B  to  J2.  (The  orientation  of  the  x  and  z  axes  are 
shown  in  Figure  2  * )  The  pressure  gradient  in  the  z  direction  decelerates 
the  negative  cross  flow  from  B  to  G  to  zero  at  some  point  between  G  and 
H.  From  H  to  Jg  the  cross  flow  is  positive  and  is  accelerated  by  the 
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* 


z  pressure  gradient.  Thus,  the  boundary  layer  fluid,  which  entered  the 
cascade  at  an  angle  larger  than  (3*,  the  free  stream  inlet  angle,  is  over¬ 
turned  and  leaves  the  blading  at  an  angle  smaller  than  p2,  the  free 
stream  exit  angle .  %.  The  process  is  illustrated  in  Figure  2. 

4.2  Limitations  in  Calculating  the  Boundary  Layer  Flow 

A  method  for  writing  the  momentum  integral  equations  in  terms  of 
four  variable®  ©x>  H,  £  and  r  is  described  in  Appendix  B.  The  shape 
factor  of  the  u  velocity  profile,  H,  was  fairly  constant  along  the 
assumed  streamline  axis,  but  the  value  of  r  varied  appreciably.. 

The  absence  of  an  expression  to  relate  r  to  the  other  variables 
prevented  any  attempt  at  a  complete  solution.  The  calculation  in 
Section  3*3  was  possible  since  the  omission  of  the  term  and 

the  assumption  of  some  constant  value  for  H  resulted  dZ  in  a  simple 
total  differential  equation  with  one  dependent  variable,  ©x. 

The  variation  of  ®x  along  the  assumed  streamline,  as  calculated  in 
Section  3-2 « 3  is  compared  with  the  measured  variation  in  Figure  l4. 


The  calculation  predicts  considerably  more. growth  in  6X  than  the 
measurements  indicate,  implying  that .neglecting  the  cross  flow  term, 

dZ 

in  Equation  (l)  is  not  justified.  This  implication  is  strengthened  by  a 
look  at  the  measured  data.  Using  the  rough  approximation  that 


i  ** 

k  " 


ze  Ze. 

(where  subscripts  refer  to  points 
in  the  flow,  Figure  6) 


we  find 


^ Z  /E 


Z  +  0.0O 7 


The  same  calculation  can  be  done  by  interpolating  the  data  for  a 
point  upstream  of  J2  (.about  0.3"  upstream  along  the  axis)  to  give 

[  ^  +  0.  0  O  G 

op 
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The  pressure  gradient  and  shear  stress  terms  at  E  and  J2  are: 


At  E- 

o ^  I  <v 
a 


—  0.008 


(measured ) 


r~~  *0  +0,0017  (Calculated  as  in  Section  5-3) 


At  J2 

'f5'  Ik  (  0,S"  (measured) 


■£.  ■hO.OOlS'  (Calculated  as  in  Section  3.3) 


The  cross  flow  term,  ~  ,  is  of  the  same  order  of  magnitude 

as  the  other  terms  contributing  to  -r~-  .  Further,  —  _  is  opposite 

o>s  oZ 


in  sign  to  the  pressure  gradient  term  and  tends  to  decrease  9X.  Thus, 
when  the  cross  flow  term  was  neglected  in  Equation  (l),  the  calculated 
growth  of  ©x  was  larger  than  that  which  actually  occurred. 


Variations  of  the  Boundary  Layer  Flow  in  the  Tangential  Direction 


The  boundary  layer  flow  along  the  streamlines,  other  them  the  one 
in  the  center  of  the  passage,  has  not  been  investigated.  The  pressure 
field  is  different  along  each  streamline  and  near  the  blades  the  com¬ 
ponent  of  flow  perpendicular  to  the  wall  no  doubt  exerts  appreciable 
influence  on  the  boundary  layer  behavior.  To  give  some  idea  of  the 
result  of  the  boundary  layer  flow  within  the  passage,  the  variation  of 
the  integral  quantities  across  the  passage  at  exit  are  shown  in  Figures  15 
and  16.  In  Figure  15,  the  large  values  of,6x*  and  9X*  near  the  suction 
side  of  the  passage  indicate  the  accumulation  of  low  energy  fluid  in 
that  region  as  a  result  of  the  cross  flow. 
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4o4  Comparison  of  the  Present  Case  to  Cascade  Tests  Made  with  a 
Uni -Directional  Inlet  Flow 

The  cascade  blades  used  in  this  study  (described  in  Section  3,2.2) 
had  been  used  in  previous  investigations  of  three-dimensional  flow  near 
the  end  walls ,  The  tests  were  run  with  an  inlet  f low,  including  the 
boundary  layer,  which  was  uniform  in  direction,  and  the  cascade  set  with; 

Stagger  angle,  §  =22.5° 

Pitch-chord  ratio,  s/c  =1.0 

The  free  stream  entered  at  Pi  =  40°  and  left  at  f3g  =  15°.  Measurements 
included; 


Velocity  profiles  at  inlet  and  exit 
Distribution  of  static  pressure  over  the  end  wall 
Distribution  of  static  pressure  around  the  blades  at  several 
distances  from  the  end  wall 


A  free  streamline  was  assumed,  and  using  the  measured  inlet  ©x,  of 
0.093  inches  and  the  static  pressure  variation  along  the  streamline,  ©x 
at  outlet  was  calculated  neglecting  the  cross  flow  term, 


bT- 

in  the  x  momentum  integral  equation,  as  in  Section  3»3«  The  calculated 
value  is  0.149  inches  while  the  measured  value  is  0 . 1^9 .  Such  agreement 
is  surprising,  for  the  cross  flow  term,  ,  at  the  exit  plane  is  of 

/u+3A  G>* 

the  same  order  of  magnitude  as >  the  x  Pressure  gradient 
term.  Apparently,  the  cross  flow  term,  tt  ,  which  is  positive  at  exit, 


changes  sign  somewhere  along  the  center  streamline  so  that  its  net  effect 
on  the  growth  of  ©x  is  small.  This  was  not  true  in  the  case  with  a 
skewed  inlet  boundary  layer  flow.  This  comparison  is  only  qualitative 
since  the  blade  stagger  angles  in  the  cascade  and  the  resulting  free 
stream  flow  differed  for  the  two  cases . 


The  distribution  of  axial  and  tangential  components  of  pressure 
force  on  the  blades  versus  y,  the  distance  from  the  end  wall  is  pre¬ 
sented  in  Figure  17,  for  the  cases  with  and  without  inlet  skewing.  The 
forces  at  the  wall  are  higher  when  skewing  is  present  due  to  the  higher 
velocities  near  the  wall  and  the  larger  angle  of j incidence  on  the  blades. 
The  forces  were  normalized  on  the  basis  of  the  free  stream. 
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5  Conclusions 

5.1  -Effects  of  Inlet  Skewing  and  Cross  Flow  in  the  Boundary  Layer 
in  the  Region.  Between  the  Blades 

The  relative  motion  between  one  blade  row  and  the  next  row  causes 
the  wall  boundary  layer  to  be  energized  and  skewed,  relative  to  the  next 
blade  row,  by  the  vector  addition  of  the  wheel  speed.  A  boundary  layer 
leaving  one  blade  row  with  large  velocity  deficiency  is  energized  and 
gets  a  "fresh  start"  as  it  enters  the  next  row,  a  very  helpful  situation 
in  a  machine  with  continuously  rising  static  pressure  along  the  flow 
path. 


The  boundary  layer  flow  at  inlet  is  initially  at  a  larger  angle 
with  the  axial  direction  than  free  stream  but  is  overturned  and  leaves 
the  blade  row  at  a  smaller  angle  than  the  free  stream.  The  cross  flow 
provides  a  mechanism  whereby  the  unusually  large  velocities  in  the 
skewed  inlet  flow  near  the  wall  can  be  diverted  and  made  more  coin¬ 
cident  with  the  free  stream  flow.  The  resulting  blade  pressure  force 
nep.r  the  wall  retains  a  fairly  high  value.  (85-90$  of  the  free  stream 
force . ) 

5.2  Designs  to  Accommodate  Wall  Boundary  Layer  Flow 

The  effects  described  in  Section  4.1  might  be  put  to  good  use.  If 
blades  were  designed  to  accept  the  skewed  inlet  flows  near  the  casing 
with  small  incidence  and  to  turn  the  flow  in  a  reasonable  rising  static 
pressure  field,  the  losses  in  the  region  should  be  reduced.  The  work 
done  by  the  blade  forces,  which  are  appreciable  (Figure  15)  would  in¬ 
crease  the  stagnation  pressure  of  the  fluid.  Blades  which  are  not  de¬ 
signed  to  accept  the  skewed  inlet  flows  are  likely  to  be  stalled  near 
the  casing.  Losbes  in  the  region  will  be  high  and  little  of  the  work 
done  on  the  fluid  will  result  in  stagnation  pressure  rise. 

5.3  Information  Required  for  Further  Analysis 

The  method  for  writing  the  momentum  integral  equations  in  the  x  and 
z  directions  in  terms  of  four  variables  presented  in  Appendix  B  must 
be  supplemented  by  additional  information  if  analytical  solutions  are 
to  be  obtained.  The  most  useful  information  would  be  two  more  ex¬ 
pressions  relating  H  and  r  to  ©x  and  £.  .  Perhaps  a  relationship 
analogous  to  the  Von  Doenhoff  and  Tetervin  (Reference  5)  expression 
for  H  in  two-dimensional  flow*  can  be  found.  Without  two  more  rela¬ 
tionships,  solutions  are  only  possible  when  H  and  r  sure  constant. 

To  know  whether  or  not  H  and  r  are  constant  for  a  particular  case 
requires  considerable  a  priori  knowledge  of  the  flow. 


=  £"'48(H'a',7S)[-f  £  f  -  -».*«.)] 
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APPENDIX  A 


Aol  Description  of  Cascade  and  Instrumentation 

The  cascade  mentioned  in  the  body  of  this  report  was  designed  to 
approximate  flow  conditions  in  the  blade  passages  of  an  axial  compressor 
rotor .  The  outstanding  features  of  this  cascade  weret 

1.  The  end -wall  boundary  layer  entering  the  plane  of  the  cascade 
was  highly  skewed;  the  flow  at  the  wall  being  at  a  high  angle 
of  Incidence. 

2.  A  moving  end  wall  ahead  of  the  cascade  was  simulated  by  use 
of  a. flow  splitter  and  a  suction  slot. 

Two  techniques  were  employed  together  to  produce  the  skewed  in¬ 
let  boundary  layer.  In  the  first  method,  a  high  velocity  jet  of  air  was 
injected  below  and  at  an  angle  with  the  main  flow.  This  jet  exerts  a 
strong  shearing  action  on  the  main  duct  flow  and  was  used  to  produce  a 
high-velocity  cross -flow  very  close  to  the  wall;  The  second  method  of 
developing  cross -flow  in  the  boundary  layer  was  to  produce  a  slight 
curvature  in  the  main  flow  upstream  of  the  cascade.  The  main  duct  flow 
was  forced  to  curve  by  suction  and  blowing  on  the  side  walls  of  the 
duct  ahead  of  the  cascade. 

Simulation  of  the  moving  wall  ahead  of  the  cascade  was  accomplished 
with  a  sharp  flow  splitter  and  a  suction  slot  in  the  end  wall  ahead  of 
the  cascade  blades.  This,  in  effect,  allows  the  velocity  adjacent  to 
the  wall  to  remain  high.  The  high  velocity  jet  passing  in  front  of  the 
cascade  is  split  by  the  sharp  leading  edge,  some  flow  goes  down  the 
suction  slot  and  some  goes  through  the  cascade  as  illustrated  in  Figure  j4. 
Figure  9  is  an  overall  photograph  of  the  cascade  taken  through  the 
transparent  side  of  the  ducting  upstream  of  the  blades.  It  shows  the 
high  velocity  jet  ducting,  the  main  flow  curvature  sucking  and  blowing 
ducts,  the  sharp  leading  edge  flow  splitter  and  suction  slot  and  the 
upstream  duct  boundary-layer -removal  slot. 

Traverses  were  made  with  a  rig  equipped  with  a  micrometer  for 
measuring  distance  from  the  wall  and  a  protractor  for  measuring  flow 
direction.  Two  probes  were  fitted  separately  into  this  rig.  The  first 
was  a  three -hole  cobra  probe  of  .014  inch  I.D.  tub'eas  which  measured 
total  pressure  and  angle .  Static  pressure  was  measured  with  an  l/8  inch 
O.D.  sphere  static  probe  set  at  the  angle  found  with  the  cobra  probe. 
Velocities  at  all  measuring  positions  were  calculated  from  the  difference 
between  measured  total  and  static  pressures.  To  measure  pressures,  probes 
were  attached  to  a  system  consisting  of  a  pressure  transducer  and  a 
galvanometer.  Pressure  distribution  on  the  blade  surfaces  were  measured 
with  the  aid  of  two  blades  equipped  with  static  taps,  one  tapped  on  the 
suction  surface  and  the  other  ontthe  pressure  side. 
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B..1  Theory  of  Turbulent  Boundary  Layer  with  Cross  Flow 

B .1.1  Momentum  Integral  Equations 

Since  expressions  relating  the  turbulent  shear  stresses  to  the  mean 
flow  are  not  known,  exact  solution  of  the  boundary  layer  equations  are 
impossible  in  the  turbulent  case.  However,  the  equations,  when  integrated 
through  the  boundary  layer,  can  be  used  for  approximate  solutions. 

Mager  (Reference  l)  has  derived  the  momentum  integral  equations  in  the 
free  streamline  coordinate  system.  When  limited  to  apply  only  to  the 
flow  along  the  x-axis  by  setting  z  =  0,  the  equations  in  a  coordinate 
system  rotating  with  constant  angular  velocity,^.,  are: 


as*  ,  i  iSr/lS, 


(2) 


(5) 


The  terms  involving are  extremely  small  for  many  cases  and  will  be 
excluded .  ' 


With  the  flow  along  the  free  streamline  just  outside  the  boundary 
layer  specified  (i  .e  ),  the  two  equations  contain  seven 

dependent  variables,  five  integral  quantities  depending  thothe  u  and  w 
velocity  distributions  and  two  wall  shear  stresses. 


The  definitions  of  the  integral  quantities  are: 


(*) 


(5) 


o 


15. 


£*r 


a  = 


d 


Obviously,  either  the  number  of  variables  must  be  reduced  or  more 
equations  must  be  obtained  to  permit  any  kind  of  a  solution. 

B.1.2  Velocity  Profiles 


It  may  be  .jpossible  to  represent  the  distribution  of  the  velocity 
components  u  and  w  through  the  b oundary* layer  by  fewer  than  five 
parameters.  The  five  quantities  ©x,  &x  ,  ©z,  6Z*  and  Qxz,  could  then 
be  replaced  by  a  smaller  number  of  variables . 

Mager  (Reference  l)  suggests  the  velocity  profiles; 

§  ■  (it 


where 


_  /  y 

V  ‘  (  s 


The  use  of  £  introduces  no  new  variables  since  it  is  the  ratio 
of  two  quantities  which  also  appear  in  the  momentum  integral  equations . 
£  seems  to  be  a  convenient  "scale  factor"  between  the  cross  flow  and 
the  main  flow,  for  _J_  w  and  u  are  of  the  same  magnitude  through  a 

6  U  U 
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good  part  of  the  boundary  layer..  If  m  said  n  are  considered  constants,  the 
quantities  defined  in  (4)  through  (8)  are  related  by: 


Constant 


(12) 


e  _L  ®*z. 

6  e,  * 


C.  ON  ST  A- NT 


-jr  =  CONSTANT 


i  =  J-  _§£-  -  CONSTA-NT 

L~  e*  a* 


(13) 


(l4) 


(15) 


and  Equations  (2)  and  (3)  can  be  expressed  in  terms  of  ©x,  6  and  TT0X . 

Mager  found  that  values  of  n  =  7,  ®  =  2  gave  velocity  profiles  which 
fitted  the  data  of  Reference  4,  but  not  the  data  from  Reference  7* 


To  determine  the  accuracy  of  assumptions  concerning  the.u  and  w 
velocity  profiles,  we  will  compare  assumed  profiles  to  measured  profiles 
from  four  sources: 

1.  Reference  4,  Point  15,  III 
Flow  through  a  90°  bend 

2.  Reference  J,  Position  la 
Flow  over  a  yawed  wing 

3.  Reference  7,  Position  lc 

Flow  over  'a  '.yawed  'wing  ■  1 

4.  Reference  6,  Figure  9>  y/s  0.5 

Wall  boundary  layer  at  the  exit  of  a  cascade  of  compressor  blades 

The  profiles  given  by  (9)  and  (10)  with  n  =  7  and  m  =  2  are  plotted 
in  Figures  18  through  21  (dashed  lines)  along  with  data  from  the  sources. 
The  assumption  of  n  =  7,  n  ■  2  is  only  accurate  for  Figure  18  and  differ¬ 
ent  values  of  m  and  n  must  be  used  to  describe  the  other  data  with  any 
accuracy.  The  profiles  given  by  (9)  and  (10),  using  values  of  m  and  n, 
which  cause  integral  quantities,  (4)  through  (8)  calculated  from  the 
assumed  profiles  to  match  the  data,  are  also  shown  in  Figures  18,  19,  20 
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and  21  (solid  lines)  along  with  the  values  of  m  and  n  usedur  The  agreement 
is  fair,  but  m  and  n  must  be  considered  variables  in  which  case  the  mo¬ 
mentum  integral  equations,  (2)  and  (3)>  contain  the  variables  Qx, 
m  and  n. 

The  use  of  5  to  non -dimens ionalize  y,  the  distance  from  the  wall, 
can  make  comparisons  between  assumed  profiles  and  measured  profiles 
uncertain,  for  5  is  sometimes  difficult  to  determine  experimentally. 

In  two-dimensional  boundary-layer  theory  (w  =  0),  the  momentum  thickness*, 
9,  which  can  be  accurately  determined,  is  often  used  to  normalize  y; 
when  cross  flows  exife$.,  9X  can  be  used.  For  the  two-dimensional  tur¬ 
bulent  boundary  layer,  Vbn  Doenhoff  and  Tetervin  (Reference  3)  suggest 
a  family  of  curves  of  u  vs.  ^  based  on  the  parameter  H,  where 

U  9 

H  -  ~  (16)  , 


The  profiles  given  by  Equations  (17)  and  (l8)  as  applied  to  data 
in  References  !<-,  6  and  7  are  plotted  in  Figures  22,  23,  2k  and  25  along 
with  the  values  of  H  and  r  used.  The  values  of  H  and  r  were  again  chosen 
so  that  the  integral  quantities  (4)  through  (8)  for  the  assumed  profiles 
matched  the  measured  values .  The  Von  Doenhoff  and  Tetervin  profiles 


of  £  vs .  ~  fit  the  data  from  the  four  sources  with  rather  good  accuracy, 

u 

while  the  agreement  in  the  cross  flow  profiles  is  fair. 

The  profiles  based  on  H  and  r  are  compared  to  measurements  made 
at  points  E,  G,  H,  I  and  J2  in  the  cascade.  Section  3*2,  in  Figures  29 
through  33.  The  agreement  is  only  fair  but  it  should  be  remembered 
that  the  accuracy  of  the  measured  data  is  poor.  The  quantity  £  is 
particularly  difficult  to  determine  since  the  accuracy  of  angle 
measurement  which  depends  on  the  dynamic  head  is  very  poor  adjacent 
to  the  wall.  The  agreement  between  the  assumed  profiles  and  those 
measured  .is  better  in  H,  I  and  J2  than  in  E  and  G.  The  cross  flow  at 
E  and  G  is  still  strongly  influenced  by  the  inlet  conditions  which 
might  be  considered  "artificial."  Further  along  in  the  flow  between 
G  and  H  the  original  cross  flow  has  been  decelerated  to  zero  and  a 
"natural"  cross  flow,  cuased  by  the  pressure  field  and  shaped  by  the 
shear  stresses  begins.  This  crude  description  of  occurrences  may 
account  for  the  differing  agreement  between  the  assumed  profiles 
(which  were  suggested  by  observations  of  the  "natural"  cross  flows) 
and  the  data. 

The  use  of  H  instead  of  n  to  describe  ~  profile,  improves  the 

accuracy  of  describing  the  data  but  complicates  mathematical  manipulation, 
for  the  analytical  relationship  between  u,  y__  and  H  is  not  known.  How- 

#  U  ©x 

ever,  the  variation  of  9xZ,  ©z  and  8Z  with  H  and  r  can  be  calculated 
numerically,  then  approximated  by  analytical  expressions.  It  Is  evident 
from  the  definitions  (4)  through  (8)  that 

i  fc*/W-aX'-  4  )"<*.)= »(".  0 
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where  g,  h  and  denote  functional  relationships. 


(19) 


(20) 


(21) 
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These  relationships  are  plotted  in  Figures  26,  27  and  28  and  it 
appears  that  g,  h  and  jj,  can  be  approximated  by  simple  functions  of  H  and  r. 

B.1.3  Shear  Stress  Lav 


In  two-dimensional  boundary  layer  calculations,  a  relationship  be¬ 
tween  'LomQf  based  on  empirical  results,  is  often  assumed.  Von  Doenhoff 
and  Tetervin  suggest  the  Squire  and  Young  formula 


(22) 


where  Re  =™ 

T 

Gruschwitz  (Reference  4)  shows  that  a  similar  friction  law 

r.  o*  o 

is  applicable  to  T  in  a  turbulent  boundary  layer  with  cross  flow  in 
a  90°  bend. 


B.1.4  Possible  Solutions 


Combining  Equations  (2)  and  (3)  with  the  shear  stress  law,  (22) 
or  (23)  and  relations  (12),  (19),  (20 ).  and  (21)  results  in  two  partial 
differential  equations  in  the  four  variables  ©x,  H,  £  and  r .  Two 
more  eqqations  relating  the  four  quantities  are  required  before  any 
solution  can  be  attempted.  The  simplest  case,  that  of  H  and  r  (or  m 
and  n)  being  constant,  has  been  treated  by  Mager  (Reference  l).  The 
profiles  in  Figure  18  were  found  representative  of  the  boundary  layer 
flow  through  a  90°  bend  and  the  relationships  (12)  through  (15 )  were 
used  with  (23)  to  express  the  momentum  integral  equations  in  two 
dependent  variables,  9X  and  &  .  For  £■  small  compared  to  1  ,  a  first 
approximation  to  the  solution  was  obtained,  while  for  £  appreciable 
compared  to  1,  a  more  complicated  second  approximation  to  the  solution 
is  presented.  The  agreement  between  calculations  and  measurements 
was  fair . 

If  H  and  r  are  not  constant,  but  their  variations  are  known, 
solutions  still  might  be  possible.  However,  an  a  priori  knowledge 
of  H  and  r  for  a  flow  is  unlikely  so  it 'seems  that  two  more  general 
relationships  between  ©x,  H,  and  r  must  be  discovered. 


Defined  by  Mager  (Reference  2). 


CONFIGURATION  AT  THE  HUB  OF  A  COMPRESSOR  WITH 
SHROUDED  STATOR  ROOTS 


FIGURE  I 
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DYNAMIC  HEAD  DISTRIBUTION  ALONG  ASSUMED  FREE  STREAM¬ 
LINE  OUTSIDE  OF  BOUNDARY  -  LAYER 


FIGURE  5 


ASSUMED  FREE  STREAMLINE 


LETTERS  INDICATE  PROBE 
POSITIONS 

NUMERALS  INDICATE  WALL 
STATIC  TAPS 

FIGURE  IS  FULL  SIZE  AND 
DIMENSIONS  MAY  BE  SCALED 
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LOCATION  OF  PROBE  POSITIONS  AND 
STATIC  WALL  TAPS 


FIGURE  6 
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BO  U  NDA  RY- LAYER  PROFILES 


FIGURE  12 
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AXIAL  AND  TANGENTIAL  FORCE  DISTRIBUTION 


FIGURE  13 


FIG-  14  -  GROWTH  OF  0X  ALONG  CENTER  STREAMLINE  OF  CASCADE  IN  REFER 
5,  MEASURED  AND  CALCULATED  (SECTION  3.3) 
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BLADE  FORCES  U.S.  SPANWISE  POSITION  WITH  AND  WITH¬ 
OUT  BOUNDARY- LAYER  SKEWING 

FIGURE  17 
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^Free  Stream 

BLADE  FORCES  U.S.  SPANWISE  POSITION  WITH  AND  WITH¬ 
OUT  BOUNDARY- LAYER  SKEWING 

FIGURE  17 


COMPARISON  OF  VELOCITY  PROFILES  ON  THE  BASIS  OF  m  AND  n  WITH  EXPERIMENTAL 
DATA  FROM  REFERENCE  4  POINT  15,  III 
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COMPARISON  OF  VELOCITY  PROFILES  ON  THE  BASIS  OF  m  AND  n  WITH  EXPERI 
MENTAL  DATA  FROM  REFERENCE  7,  POSITION  la 


DATA  FROM  REFERENCE  6  FIGURE  9 


COMPARISON  OF  VELOCITY  PROFILES  ON  THE  BASIS  OF  H  AND  r  WITH  EXPERI 
MENTAL  DATA  FROM  REFERENCE  6,  FIGURE  9 
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DATA  FROM  REFERENCE  7  POSITION  Ic 


FIG.  29  COMPARISON  OF  VELOCITY  PROFILES  ON  THE  BASIS  OF  H  AND  r  WITH 
EXPERIMENTAL  DATA  FROM  REFERENCE  5. 


